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Using a new scheme of the derivation of the non-hnear a- model we consider the electron motion in 
a random magnetic field (RMF) in two dimensions. The derivation is based on writing quasiclassical 
equations and representing their solutions in terms of a functional integral over supermatrices Q 
with the constraint — 1. Contrary to the standard scheme, neither singling out slow modes nor 
saddle-point approximations are used. The cr-model obtained is applicable at the length scale down 
to the electron wavelength. We show that this model differs from the model with a random potential 
(RP). However, after averaging over fluctuations in the Lyapunov region the standard a- model is 
obtained leading to the conventional localization behavior. 
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I. INTRODUCTION 

Description of the two dimensional {2D) electron mo- 
tion in a random magnetic field (RMF) is of a consid- 
erable interest for both experimentalists and theoreti- 
cians. Two dimensional electron systems in a random 
magnetic field were realized in a number of recent exper- 
iments when a high-mobility heterostructure was located 
under an overlayer with randomly pinned flux vortices in 
a type-II superconducting gate^ or type-I superconduct- 
ing grains^ or a demagnetized ferromagnet'^. From the 
theoretical point of view the RMF model is an example 
of a system with the interaction which is realized through 
an effective gauge field. In particular, this model arises in 
the theory of quantum Hall effect with a half-filled Lan- 
dau level"*. Another application of this model is a gauge 
field description of the doped Mott insulators^. 

One of the most important problems in the RMF mod- 
els is the question about localization of electron states. 
This question has been studied in many numerical works 
and very different conclusions were drawn: from a) all 
the states are localized, Refs.^~^ to b) there may be a 
band of delocalized states Refs.^^*** and c) all the states 
are localized except those with the precisely zero energy, 
j^gfgi5,i6 rpj^g problem of comparison of the results 
obtained in different numerical calculations is a quite 
complicated task partly because extended states and the 
states with very large localization length can very often 
be hardly distinguished from each other. 

From the point of view of the generally accepted scal- 
ing theory of localization*^ the RMF model should not 
be different from the model describing the electron mo- 
tion in a random potential in a homogeneous magnetic 
field. In both the cases all electron states are expected 
to be localized in 2D in an arbitrarily weak random po- 
tential. Using the supersymmetry technique*^ this pre- 
diction was checked in several works by deriving a proper 
CT- model. The authors of Ref.*^ used the standard scheme 



of the derivation finding first the saddle-point in the in- 
tegral over supermatrices Q and expanding then in slow 
modes near this saddle point. As a result, they obtained 
a standard diffusive unitary a-model similar to what one 
has for the model with a random potential (RP) and the 
broken time reversal symmetry. The long range charac- 
ter of correlations of the random vector potential, which 
is possible even if the correlations of the magnetic field 
are short ranged, did not play any role. 

A possibility of a new term in the cr-model due to spe- 
cial character of the correlations of the vector potential 
was discussed later in Refs.^""^^. This was done by con- 
sidering more carefully short distances. A ballistic a- 
model similar to that of Ref.^^ was derived in Refs.^"'^^ 
and the calculations were checked by direct diagrammatic 
and path integrals methods^*. The final conclusion of 
these works was that the cr-model maintained the stan- 
dard form*^ corresponding to the unitary ensemble unless 
the correlations of the magnetic field were long ranged. 
This was considered, as usual, as the proof of the localiza- 
tion. An additional term in the a-model was still possible 
if the correlation of the magnetic field was proportional 
to Ref.^^, where q is the momentum, and this could 
lead to antilocalization (see also Ref.^**). However, no 
possibility to obtain anything but the standard unitary 
CT-model and, hence, the localization for any finite range 
correlations of the magnetic field was seen finally from 
these works and no difference between the RMF model 
and the RP model with a magnetic field was found even 
in the ballistic case. 

Nevertheless, the question about the localization in the 
RMF model in 2D was raised again in a recent numerical 
work^''. On the basis of the numerical study the author 
of Ref.'^^ suggested quite a different scenario of the elec- 
tron motion in the RMF model arguing that there could 
be some "hidden degrees of freedom" that lead to essen- 
tial deviations from the standard scaling description of 
disordered systems. 
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This resuh chaUenges the analytical results obtained 
on the basis of the cr-model description but it is fair to say 
that the previous analytical study was not complete. All 
calculations were carried out using the traditional form of 
the ballistic cr-model^^'^^'^^ with a conventional collision 
term. However, this form may be used for a long range 
disorder at sufficiently long distances only. The deriva- 
tion of such a cT-model is based on finding a saddle point 
in the integral over the supermatrices Q and expanding 
in slow modes. This procedure fails at short (but still 
much exceeding the wave length \p) distances. As a re- 
sult, the form of ballistic cr-model is not applicable at the 
lengths smaller than a characteristic length ^ Xp and 
this puts doubts on some conclusions drawn previously. 

The saddle-point approximation is equivalent to the 
self-consistent Born approximation (SCBA) and cannot 
be good for a long range disorder. At the same time, 
even short range correlations of the magnetic field corre- 
spond to long range correlations of the vector potential 
and this problem is inevitably encountered in the RMF 
model. The diagrammatic expansion of Ref.^^ also starts 
with the SCBA for one-particle Green functions and one 
encounters the same problem. 

In order to circumvent the problem related to the 
use of the saddle point approximation and the expan- 
sion in the slow modes we suggested recently another 
scheme^''. This method is based on equations for quasi- 
classical Green functions and resembles the phenomeno- 
logical approach of Ref.^'^. However, in contrast to the 
latter, we do not average over disorder in the begimiing of 
the calculations and do not decouple an effective interac- 
tion by integration over an auxiliary field. Our approach 
is exact in the quasiclassical limit and a resulting ballistic 
tr-model is applicable at all distances exceeding the wave 
length Ai? . It can be reduced to the conventional ballistic 
(T-model after a coarse graining procedure and the latter 
is applicable at distances exceeding a Lyapunov length 
II introduced in Ref.^®. At distances smaller than II the 
form of the term due to disorder is different from the 
standard collision term. 

In Ref.^^ we derived the ballistic cr-model for the RP 
models and now we present an analogous derivation for 
the RMF models. It turns out that the terms in the 
ballistic cr-models describing the disorder in the RP and 
RMF models differ from each other. They can become 
similar only after carrying out the coarse graining proce- 
dure. We show that this procedure can be performed in 
the same way as for the RP problem, which leads to a 
similar reduced cr-model. 

The paper is organized as follows: In the Chapter II 
we introduce a partition function generating correlation 
functions of interest in terms of a functional integral over 
supervectors ijj. We derive equations for Green function 
and simplify them using a quasiclassical approximation. 
Introducing quasiclassical Green functions we rewrite the 
equations in a gauge invariant form. The solution of the 
equations is found in terms of an integral over superma- 
trices Q with the constraint = 1, which allows us to 



average over the RMF. 

In the Chapter HI we integrate over fluctuations in the 
Lyapunov region and come to a reduced cr- model with a 
collision term. 

In the Appendix we consider the problem of the cor- 
relation of two particles moving in a RMF and find the 
characteristic time of this correlation. 



II. FORMULATION OF THE PROBLEM. 
QUASICLASSICAL APPROXIMATION AND 
DERIVATION OF THE <r-MODEL. 

In the present work we follow the method of deriva- 
tion of the cr-model suggested in our previous work^''. In 
order to make the presentation self-contained we repeat 
the main steps of the derivation. 

We start our consideration with the introduction of the 
partition function Z[a\ 



Z[a] 



J exp(- 



■LM)Di; 



(2.1) 



LabP] = -i I ip{r) ( H{r) - £ + | + ^^^^^A ) V'(r)cir + 



+i J il){r)a{r)il){r)dr 



where ip are 8-component supervectors^^ and the Hamil- 
tonian H{r) in Eq.(2.1) is taken in the form 

^(r) = (-iVr - ^^3A(r)) ' /2m - + ^(r) (2.2) 

The last term in Eq.(2.1) contains a source function 
a(r). Choosing this function in a proper form and taking 
derivative in it one can obtain correlation functions. For 
example, the level-level correlation function R{lo) can be 
written as: 

1 1 

R{oj) = - - lim Re— — — Z[o] (2.3) 

^ ^ 2 2(7rz/F)2 ai=a2=o daida2 

where the source a(r) is the following matrix: 



a(r) = 









-Q!2 



a,,2 = ^{l-k) (2.4) 



Here k is the diagonal matrix with elements ±1 in 
fermionic and bosonic blocks respectively^^. 

The Hamiltonian H{r), Eq.(2.2), contains both scalar 
and vector potentials w(r), A(r) that are assumed to be 
random functions of the space coordinates distributed 
according to the Gauss law, fa is the third Pauli matrix 
in the particle-hole space. Below we consider a general 
case when the scalar potential u(r) contains both the 
short range Us (r) and long range ui (r) parts with the 
characteristic correlation lengths of the order and larger 
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than the Fermi wavelength Xp = {2TrpF)~^ respectively. 
Their statistics are determined by the pair correlation 
functions: 



M(r) = CM^{r)C^ 



{us{r)us{r')) 



1 



{ui{r)ui{r')) ^W{r~r') 



S{r - r') 



(2.5) 
(2.6) 



where the function M^(r — r') is assumed to fall off over 
a length Xp- Statistics of the magnetic field will be 
introduced later. Although the main goal of this paper is 
to study the RMF model, we add the scalar potential into 
the Hamiltonian for a more explicit comparison between 
the RMF and RP models. 

Following the standard approach of Ref.^^ one would 
average the partition function Z[a], Eq.(2.1), over the 
random external fields and then, singling out fluctuations 
slowly varying in space and integrating over an auxiliary 
smooth matrix field Q, decouple the interaction term 
(tp'tp)'^ that appears after the averaging. This method 
was recently used, e.g., in Ref.^^ in a derivation of the 
ballistic a-model for quantum billiards and in Ref.^°'^^, 
where the two-dimensional electron gas was considered 
in a random magnetic field. As it has been mentioned in 
the section I the latter problem is rather specific because 
the vector potential A(r) can have long range correla- 
tions even if correlations of the magnetic field are short 
ranged. 

The singling out of slow modes with the subsequent 
decoupling of the interaction by integrating over an aux- 
iliary smooth matrix Q is not a rigorous procedure be- 
cause some part of the interaction is assumed to be irrel- 
evant and is neglected. Although this assumption works 
well for short range impurities, it is not justified for long 
range correlations. Below we use another method based 
on the Green function and quasiclassical approximation 
of Ref.^^. This method allows one to derive a cr-model 
applicable down to the length scale of the order of the 
wavelength Xp. 

Following Ref.^'' we average over the short range po- 
tential Us{r), decouple the interaction term appearing 
after this averaging using the standard integration over 
an auxiliary smooth matrix field M(r) and finally rewrite 
the partition function as follows: 

Z[a\ = j Zi[J]exp^-^ j StvM'^{Y)d^DM (2.7) 

where 

Zi[J] = j exp(-Lj[^])£»V (2.8) 

The Lagrangian coincides with Lalip], Eq.(2.1), 

provided the substitutions Ug{Y) = and M(r) J{r) = 
ia{r) + M{t)/2ts are made in the Lagrangian La[il}\, 
Eq.(2.1). The structure of the matrix M(r) can be 
foimd in the book, Ref."'^*. It is important that M{y) 
is self-conjugate: M(r) = M(r) where the bar means the 
"charge conjugation" 



C = A( 



ci 

C2 



Cl 



-1 

1 



C2 



1 

1 



(see also Ref.^^). 

Following Refs.^^'^^ we introduce the Green function 
G(r,r') 

Gfa/3(r,r') = Zf i[J] j V"(r)i/;'5(r')e-^^[^]£'V (2-9) 

For the most correlation functions of interest the source 
function a(r) can be chosen to be self-conjugate. If this 
is the case the Green function satisfies the equation 



H{v)-e+- + ^^K + iJ{v) 



G(r,r') =i(5(r-r') 
(2.10) 



Eq.(2.10) was previously studied in the absence of the 

magnetic field in the quasic;lassical approximation using 
a method of a quasiclassical Green function, Refs.^^'^''. 
This method is based on the assumption that the exter- 
nal fields and sources are smooth functions (i.e. slowly 
changing over the wavelength Xp). Within this method 
the Green function G(p,R) can be rewritten using the 
Wigner transformation 



G(r,r') 



s*p('-'"')G(p,R), R=(r-Fr')/2 



The function G (p, R) has a sharp peak at the Fermi sur- 
face p = pfh. This property is due to the fact that the 
long range fields and sources weakly disturb the shape 
of the Fermi surface. Integrating the Green function 
G(p, R) over the absolute value of the momentum p re- 
sults in a new function gn{r) that depends on the centre 
of mass coordinate R and the unit vector n ~ p/p deter- 
mining the direction at the Fermi surface. The coordi- 
nate dependence of this function turns out to be smooth 
and therefore .gn(r) may be considered as the quasiclas- 
sical approximation of the exact Green function G(r, r'). 
On the other hand, the partition function Zi [J] , Eq. (2.8) , 
can be expressed through (?n(r). 

Before we start the calculation following this procedure 
let us make some remarks about differences between the 
RP and RMF models. First, the presence of the magnetic 
field breaks the time-reversal symmetry and, hence, ex- 
citations sensitive to the time reversal arc suppressed. 
Therefore we consider only such correlation functions 
that can be obtained from the sources a(r) commuting 
with T3. The part of the Green function anticommuting 
with T3 is negligible and may be omitted from the further 
consideration. 
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The second remark is related to the physical aspects 
of the quasiclassical approximation in the presence of 
a magnetic field. It is known that systems placed in a 
magnetic field are invariant with respect to the magnetic 
translations = exp [(Vr — z(e/c)f3A) a] instead of the 
ordinary ones'' ^. The difference between these transla- 
tions is relevant for an infinite system even if the mag- 
netic field is weak. This means that electron states are 
to be characterized not by the ordinary momentum pkin 
determining the kinetic energy but rather by the gener- 
alized momentum p = p^m + (e/c)f3A(r). The general- 
ized momentum p is a well-defined quantum number if 
the magnetic field is weak: 



vf 

th > Af, th = — 



(2.11) 



where coh = eH/mc is the Larmor frequency and vp- 
Fermi velocity. Inequality (2.11) coincides with the con- 
dition of the applicability of the quasiclassical approxi- 
mation. The Fermi surface is defined in the space of the 
generalized momentum p and, contrary to the case of 
zero magnetic field, has a rather complicated form. The 
value of the momentum p at the Fermi surface strongly 
depends on the direction n = p/p. Therefore we change 
the definition of the quasiclassical Green function by re- 
placing the integration over the absolute value of the gen- 
eralized momentum p by that of the kinetic one pfei„ (see 
e.g. Refs.30): 

5n(r) = ^jdiG(p+ -hA{T), r) (2.12) 

where the function G(p, r) is the Green function taken 
in the Wigner representation and ^ = p^/(2m) — £f, 
n = p/p. The quasiclassical Green function gii{v) de- 
fined by Eq.(2.12) is gauge- invariant. The logarithmic 
derivative of the partition function Zi[J], Eq.(2.8), can 
be estimated as follows: 



6\nZi[J] 



^G(r,r) 



^ j 9n{r)dn 



(2.13) 



where is the density of states at the Fermi surface. Per- 
forming the Wigner transformation we subtract Eq.(2.10) 

from the conjTigatcd one, then integrate the result over 
^ as in the Eq.(2.12) and obtain in the quasiclassical ap- 
proximation: 



[^UFnVr + — f3B{r)d^ - pp Vru(r)a„) gn{r)- 



+ id) 



[A, 5„] - [J(r), Qn] = 



(2.14) 



In Eq. (2.14), B{r) = d^Ay — dyA^ is the magnetic field, 
5n = ^ipd,^ TO., = (—sin 95, cosip). In this approxi- 
mation the solution of the Eq.(2.14) is to be sought with 
the usual constraints^ 



9l{^) = 1 
and the boundary condition 



es 



(2.15) 



(2.16) 



where r € S stands for points on the surface of the sam- 
ple and n_L means the component of the vector n per- 
pendicular to the surface. Following Ref.^^ we write the 
solution of Eq.(2.14) in terms of a functional integral over 
supermatrices Qn (r) 

<?„(r) = Z^'[J] f g„(r)exp (-^$^[Q„]) DQ„ 
Z2[J]= [ exp(-^a>j[Q„])DQ„ (2.17) 



*j[Qn] = Str / drdn 



AT„(r)(t;FnVr + ^^f35^ 



mc 



p^iVru(r)V„ )r„(r) + ( !l^^i±l^A - J(r) ) Q„(r) 



(2.18) 



Q„(r) = r„(r)AT„(r), f'„(r)T„(r) = 1 

where stands for the derivative in the angle. The inte- 
gration in Eq.(2.18) is performed over the self-conjugate 
supermatrices 

<3n(r) = Q„(r), Q„(r) = CQ^„(r)C^ 

with the constraint Q\(j) = 1 and 



<3n^(r)|s=Q-n^(r)|s 



(2.19) 



at the surface S of the sample. The structure of the 
supermatrix Qn coincides with the structure of the su- 
permatrix M(r). We do not demonstrate here the equiv- 
alence of the matrices g„ (r), Eqs.(2.12), (2.18), and re- 
fer to the proof given in Ref.^^. We mention here only 
that both the matrices are the logarithmic derivatives in 
the matrix J(r) of the partition functions Z\\J\, Z2[J] 
respectively. Hence, these functions are equal to each 
other up to some factor that is independent of J(r). Due 
to the supersymmetry Zi[J] = Zi\J\ = 1 for J(r) = 0, 
which means that the factor is unity and the partition 
functions are equal to each other 



Zy{J\=Z2\A 



(2.20) 



Below the magnetic field -B(r) is considered as a ran- 
dom function with a Gaussian distribution and the pair 
correlation function of the form 

(B(r)B(r')) =2(^)'a;2W^s(r-r') (2.21) 
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where is lUc is a coefficient that has a meaning of the 
characteristic frequency of the cyclotron motion and the 
function Wb{i^ — r') is assumed to fah off at distances 
|r — r'l > b and to be normalized as Wsir = 0) = 1. 
The length b characterizes the decay of the correlations 
of the RMF B (r). Substituting Eq.(2.20) into Eq.(2.7) 
and averaging the result over the magnetic field and long- 
ranged potential u;(r) we find for the partition function 
Z[a] Eq.(2.7) 

Z[a] = J eM-F[Qn])DQn (2.22) 

where the free energy functional F[Qa\ has the form: 

F[Qn\ = FkiniQr.] + -Fimp[Qn] + F^^M + F„[Q„] 



Nevertheless, at longer distances the RP and the RMF 
models arc very similar and we show this in the next 
Chapter carrying out a coarse graining procedure sug- 
gested in Ref.^^. The latter means integrating out de- 
grees of freedom at distances inside the Lyapunov region. 

For simplicity of the presentation we will consider in 
the next Chapters only effects related to the random 
magnetic field and disregard the scattering on the ran- 
dom potentials omitting F^rnp [Qn], -Fimp[Qn] in the free 
energy Eq.(2.23). Accordingly, we will consider the sym- 
metry of the supermatrices Q corresponding to the uni- 
tary ensemble. We will study the behavior of the a- 
model, Eq.(2.23), on different length scales and discuss 
the connection of this model with the models previously 
obtained in Refs.^^~^^. 



FkinlQa] = — Str / drdn 



Ar„(r)vj^nVrr„(r) 



+,f5^A_a)Q„(r) 



^^mp [On] = -l (^^) I drdndr'dn'\7lVi,W{r - r') 
xStr[Ar„ (r) V^,r„(r)]Str[Ar„,(r')V^.,r„,(r')] (2.23) 



FmiQn] =[y''' 



drdndr'dn'WB{r-r') 
Str (Af3T„(r)i5^T„(r)) Str (Af3T„,(r')ia^'T„,(r')) 



and 



V„ 



(2.24) 



The first term Fkin [Qn] describes the free motion and 
is what remains when external fields and impurities are 
absent. The second and the third terms FimplQn], 

^implQ^] ^'"^ responsible for the scattering on the long- 
and short-ranged potentials respectively. The last term 
Fm [Qn] is due to the presence of the random magnetic 
field. Correlation functions of interest can be obtained by 
calculating derivatives in the source d(r) of the partition 
function Z[d], Eq.(2.22). 

It is important to emphasize that the structure of the 
terms -Fimp[Qn] and F,m[Qii] describing the electron scat- 
tering on the random potential and on the random mag- 
netic field, respectively, is clearly different. The term 
Fimp[Qn] contains the components of the gradients par- 
allel to the plane, whereas the term Fm[Qn] contains the 
perpendicular one. 



III. REDUCED fT-MODEL 

The (T-model obtained in Eq.(2.23) is valid for the 
length scales down to the wavelength Xp and has the form 
which differs from the a- model found in the Ref.^°'^^. 
The latter model has been derived for the spatially un- 
correlated magnetic field and is applicable at the length 
scale restricted from below by the single-particle relax- 
ation length I but not by the wavelength A_f. The length 
I could not be consistently estimated within the consid- 
eration of Refs.^"'^^ and remained without a clear phys- 
ical interpretation. At the same time, the analysis of 
Refs.^®'"^^'^^ leads to the conclusion that the role of this 
length is played by the Lyapunov length II = vpT^. Here 
tl is the inverse Lyapunov exponent and is the time dur- 
ing which two close trajectories increase the distance be- 
tween them by a factor of the order of unity. On the 
other hand, according to the Ref.^^, Tp is the time which 
is required for two scattered particles to diverge over the 
distance of the order of the range of the potential (or 
the correlation length). In the Appendix we discuss the 
problem of the particle motion in a RMF and estimate 
the Lyapunov length II for weak fields as 



II 



Hr 



2/3 



(3.1) 



This result shows that the Lyapunov length / l is between 
the correlation b and transport Itr lengths: b <^ II hr- 
The Lyapunov length 1^ divides the length scales into 
two regions. At small distances, two particles propagate 
in the same magnetic field and correlations between them 
are relevant. Following the terminology of Ref.^* we call 
these distances the Lyapunov region. In the second re- 
gion when the scales of interest are larger than the Lya- 
punov length, the motion of the particles is not correlated 
and they are scattered by the RMF independently. This 
can be called the collision region because the correspond- 
ing classical motion at such distances is described by the 
conventional Boltzmann equation with a collision term 



5 



corresponding to the scattering on the RMF. The elec- 
tron motion at these long distances should be described 
by a reduced cr-model and one can expect that this re- 
duced cT-model is just the a-model of Ref.^"'^^. In order 
to obtain the reduced cr-model one should integrate out 
in Eqs. (2.22, 2.23) the degrees of freedom related to the 
Lyapunov region. This coarse graining procedure has 
been worked out in Ref.^^ for the RP model and we will 
repeat it now for the RMF model. 

First, one should explicitly decouple the original mode 
Tn(r) into the "slow" and "fast" parts. We make this sep- 
aration in the way preserving the rotational invariance of 
the initial model Eq.(2.23): 

T„(r) =r„(r)y„(r) (3.2) 

Here Tii(r), Vn{r) are "slow" and "fast" modes describing 
the fluctuations in the collision and Lyapunov regions re- 
spectively. As soon as the mode separation is made one 
should substitute Eq.(3.2) into the free energy F[Qn], 
Eq.(2.23), and then average it over the "fast" fluctua- 
tions Vn{r): 

Z[a] = [ e-^-ff^'^-^DQn (3.3) 

where 

^-F.fAQ.] = J eM-F[Q'^^]-Fint[Q^^\Qn])DV^ 

(3.4) 

QW(r) = F„(r)Ay„(r), g„(r) = T„(r)A^„(r) 

The functional i^[(5n in Eq.(3.4) coincides with the free 
energy Eq.(2.23) provided the source is omitted in the lat- 
ter expression. The functional Fint [Qn'' , Qn] determines 
the interaction between the fast and slow modes Qn ' , Qn 
and has the form: 

F'int [Qii ) Qn] = -Pfein [Qn ^ J Qn] + Fm [Qii ^ i Qn] 



^n(r) = Tn(r)ia^T„(r) 

Before the averaging over the fast fluctuations Qn ^ we 
make the following essential remark. 

The separation into the fast and slow modes, Eq.(3.2), 
requires a more acciirate definition. The point is that 
the excitations in the model Eq.(2.23) reveal a strong 
anisotropy in the phase space (r, n) due to the specific 
form of the free energy functional, Eq.(2.23). Since only 
the first order derivatives in r and n enter the free energy. 
Eq.(2.23), the dependence of the excitations on the co- 
ordinates (r, n) will resemble a propagation along a clas- 
sical trajectory. Such an anisotropy demands a care and 
should be performed in an invariant way. As in Ref.^^, 
the scale separation can be performed introducing an ad- 
ditional term into the functional i^[Qn^], Eq.(3.4), 

FL[Q^n^] = -^AiStr j drdnAQ^°\r) (3.6) 

Then, we extend the region of the integration over Qn ' (r) 
to all possible matrices with the constraints Eq.(2.19). 
The parameter Al is just the Lyapunov exponent t^^ 

and the term F^lQn^], Eq.(3.6), serves to suppress fluc- 
tuations of the matrices Qn^ outside the Lyapunov re- 
gion. 

As soon as the mode separation is properly defined one 
can carry out the integration in Eq.(3.4) and evaluate the 
effective energy -F'e//[Qn]- We perform this computation 
using the cumulant expansion in Fint, Eq.(3.4) and ap- 
proximation of the weak magnetic field. In the same way 
as it was done in Ref.^'' for the model of the long-ranged 
disorder one can show that this is an expansion in powers 
of the operator ^^Vr which is small outside the Lyapunov 
region. Considering only the first order we find 

Feff[Qn] = {Fint[Q^^\Qn])o (3.7) 

where the brackets (. . .)o stand for integration over 
Due to the supersymmetry {Qn\'r))o = A, which gives 



nin{Q^n\Qn] = ^StT / dvdn 



2 
+i 



{FLn[Q'^\Qr.])0 = Fkin[Qn] 



u) + i5 



a(r)')fn(r)QL°^(r)Tn(r) 



FmiQn J Qn] — 



drdndr'dn'Wsir -r') 
X Str [fsQi") (r)$n(r)] Str [tsQ^^) (r')$n' (r')] + 
+2 (Jy^^T J 'i^dndr'dn'WBir - r') (3.5) 
X Str [f3QL°) (r)$n(r)] Str [AfaFn' {r')id^'Vn' (r')] 



(3.8) 



with the same functional Fkin[Qn] as in Eq.(2.23). The 
second term in the functional -Fm[Qn"\ Qn] Eq.(3.5) van- 
ishes after the averaging due to the symmetry as well. 
The contribution coming from the first term can be di- 
vided into two parts: the first one comes from the re- 
ducible average and coincides with the magnetic energy 
Fm[Qn] of the initial functional, Eq.(2.23), whereas the 
other is given by the irreducible average {{Qn^Q^}))o = 
{Q'^^Q^n'h ~ (Qi°'>o(Qi°')o of the supermatrices Q^°\ 

In order to find the contribution coming from the irre- 
ducible average we consider the matrix 
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5ni(ri;a) 



(QL?(ri)exp 


^Str/rfrdna„(r)QL°^(r) 


)o 


(exp 


'^StrJdrdnan{r)Q^n\r) 


>o 



(3.9) 



where the new source dn(r) is 

an(r) = a;(r)f3$„(r), 

a(jc) is some function. Due to the supcrsymnictry 
g„(r; a = 0) = A. The first derivative in the function 
a(r) gives 



Sa{r2) 



a(r)=0 



(3.10) 

xStr J dn'Q^°}{r2)h<^n'ir2)))o 



On the other hand, the matrix gn{r; a) satisfies the equa- 
tion 

VFn\/r9n{r;a) +t [A, gn{r;a)\ = 



Q!(r)[f3$„(r), ffn(r;a)] 



(3.11) 



and condition g^{r;a) = 1. Differentiating in a(r) both 
sides of this condition and then putting a{r) =0 we find 
that the matrix Sgn{r;a)/Sa{r')\a=o in Eq.(3.10) is off- 
diagonal. Eq(3.11) can be considered for the off-diagonal 
part of the matrix gn{r;a) and rewritten in the integral 
form 

g^{r;a)= j dr'g„(r - r')Q(r') [f3$„(r), ff„(r;Q)]^ 

(3.12) 

where the superscript ± stands for the part of the super- 
matrices anticommuting with A. The kernel ^n(r — r') is 

the solution of the equation 

[vpr^Vr + i(w + iAi)A]e„(r - r') = 5{r - r') (3.13) 

Differentiating in Q:(r) both sides of Eq. (3.12) and 
putting a(r) =0 we obtain 



~2 



((QW(r,)Str j dn'Qi°;(r2)f3$„.(r2 



)»o 



ani(ri-r2)[f3$„,(r2), A] (3.14) 
Substitution of Eq. (3.14) into the Eq.(3.5) gives 

{F'^[Qr., QL°^])o = FmlQr.] - T^vu?, J dvdv' dnWeiv - v') 
xStr [$^(r)g„(r - r')A$^(r')] (3.15) 



Characteristic values of the difference r — r' in ^n(r — r') 
are in the Lyapunov region, whereas 'l'n(r) is a smooth 
function. This allows us to make the replacement r' ^ r 
in one of the in Eq.(3.15). The integral over the differ- 
ence /9 = r — r' is calculated as follows. First, we rewrite 
this integral using integration in the momentum space 
instead of the coordinate one 



/ 



Gr.{p)WB{p)dp 



(2^ 



WB(q) 



(3.16) 



The momentum q may be considered as the transfer mo- 
mentum q = p' — p, where p' = ppn'. p = ppn are 
momenta of a particle after and before the scattering. 
Since for a weak scattering the characteristic length b of 
the distribution VFB(r — r') is much smaller than the Lya- 
punov length, b <^Il, Eq.(3.1), the fraction in Eq.(3.16) 
can be replaced by the ^-function 



Uirnq — {u + iAi,)A 



dq 



Wb (q)^(t;Fnq) 



(3.17) 



The ^-function fixes the value of the final momentum 
p' on the Fermi surface: (5(wi?nq) = S[vpn{p' — p)] = 
d[{ds/dp){p' — p)] = (5[£(p') — £(p)]. Integrating over 
the energy e' = e{p') we find for the integral, Eq.(3.16), 
the following expression 

-TTJ/A J dn'WB\pFin - n')] (3.18) 

Taking together Eqs.(3.15), (3.16), and (3.18) we obtain 
the free energy F^f / [Qn] of the reduced ct- model 



Feff[Qr.]=F[Qn]+F'[Q^], 



(3.19) 



FlQr.] 



drdnStr 



fn{r)vFnVrTn{r) + (3.20) 



T 



z I ^^A - a{r)\ Q„(r) + -^^{dM 



F'[Qn] = - (^t^c) ^ J drdndr'dn'WB{r - r') (3.21) 

xStr (Af3T„(r)a^T„(r)) Str {Ahfn'ir')d^'Tn'{r')) 

The collision term in the free energy functional is ex- 
pressed through the transport time Ttr 



2TTVTtr) ^= j du'LO^WBlppin' -n)] (3.22) 



and agrees with the results of the Refs.^"'^^'^^ where the 
RMF and long-range disorder models, respectively, were 
considered in the limit of small scattering angles. The 
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second term J^'[Qn] in Eq. (3.19) is small and can be ne- 
glected. This can be easily understood using the fact that 
the Fourier transform of the function Wb in Eq. (3.22) 
contains momenta of the order oi pp, which corresponds 
to short distances of the order of Xp- In contrast, the 
main contribution to the integral over the coordinates in 
Eq. (3.21) comes at weak RMF from larger distances of 
order II where the function Wb is small. Therefore, ev- 
erywhere below wc will imply that the reduced ballistic 
cr- model is described by the free energy functional -F[(3n] 
fromEq. (3.20). 

Thus, wc have demonstrated that, although the bal- 
listic cr-modcl for the RMF is different from the one for 
the RP (the terms F,„p[Q„] and F,n[Qn] in Eq. (2.23) 
are different), the reduced cr- models describing the elec- 
tron motion exceeding the Lyapunov length II have the 
same form of Eq. (3.20). The similarity of the RMF and 
RP models has been emphasized in Ref.^^ and the final 
conclusion of Ref.^^ was the same. However, the meth- 
ods used in these works were based on writing first the 
self-consistent Born approximation for one particle Green 
functions (saddle point equation in the cr-model formu- 
lation) and on a subsequent expansion in slow modes, 
which could not be justified at short distances. Now 
we see that the equivalence of the RMF and RP mod- 
els can hold at distances exceeding the Lyapunov length. 
This naturally leads to the equivalence of the diffusive 
cr-models that can be written in the standard form 



F[g] = ^Str 



D{VQ) +2i{uj + iS)AQ dr (3.23) 



where D = VpTtr/2. For the RMF problem the transport 
time Ttr is given by Eq. (3.22). 

Eq. (3.23) is valid unless the correlations of the mag- 
netic field are very long ranged. Only if 



{B,B_,) ~ q- 



(3.24) 



an additional term can appear^^. The symmetry of the 
difi'usive cr-model, Eq. (3.23), corresponds to the imi- 
tary ensemble and one comes to the standard conclusion 
about the localization. 

Of course, the coarse graining procedure leading to 
the ballistic a-model, Eq. (3.20), is possible only if 
the ground state of the initial cr-model, Eq. (2.23), is 
achieved at Q = A. One can imagine such functions 
Wb (r — r') that this ground state is no longer stable. 
However, this could be possible only if the Fourier trans- 
form Wb (q) was negative for certain q, which is excluded 
in the case of real magnetic fields. Therefore, beyond the 
Lyapunov region, the ballistic cr-model, Eq. (3.20), and. 
correspondingly, the diffusive a-model, Eq. (3.23), seem 
to be unavoidable. 



IV. DISCUSSION 

In the present paper we considered the problem of 
the two-dimensional electron gas in a random magnetic 



field (RMF) using the non-linear supermatrix cr-model 
approach. We derived a ballistic cr-model avoiding the 
standard scheme based on finding a saddle point in the 
integral over supervectors and expanding in slow modes 
near this point. Such a scheme explicitly relies on the 
assumption of a sufficiently short correlation length of a 
random potential (see e.g. in Ref.^^) and its validity for 
a long range disorder is not clear. As the vector potential 
entering the RMF model has a large correlation length 
even when the magnetic field is (5-correlated in space, the 
procedure of singling out slow modes used in the stan- 
dard derivation is not well justified at least at not very 
large distances. Besides, the saddle-point approximation 
is hardly allowed in this case as well. 

Instead of following the standard scheme wc used the 
method based on writing quasiclassical equations for 
Green functions and the exact representation of their so- 
lutions in terms of integrals over supcrmatrices with 
the constraint = 1. This method needs neither sin- 
gling out the "fast" and "slow" parts from the interac- 
tion nor the saddle-point approximation. Conditions of 
the applicability of the method coincide with those of the 
quasiclassical approximation. Therefore, the cr-model ob- 
tained should be applicable over the distances down to 
the Fermi wavelength, which makes it more general in 
comparison with the cr-models derived earlier on the basis 
of the standard scheme, Refs.^^'^^'^"'^^. The latter mod- 
els are justified at distances exceeding the single-particle 
mean free path I as in the Ref.^°'^^ or the transport length 
Itr as in Ref^^. 

We have demonstrated that similar to the problem of 
long range random potential, there is a characteristic or 
Lyapunov length dividing the length scale into the 
Lyapunov and collision regions. The first region corre- 
sponds to the small distances over which the particle mo- 
tion is strongly correlated. Correlations disappear over 
the larger lengths where the particle interaction can be 
considered in terms of collisions. In the Appendix we es- 
timate the Lyapunov length for RMF problem restricting 
our consideration by the limit of a weak field. The esti- 
mated length is expressed through the transport length 
Itr and the correlation length b of the RMF by a formula 
similar to the one obtained previously in Ref.^* in the 
model of a long ranged potential. 

The reduced cr-model obtained after the integrating 

over the fluctuations in the Lyapunov region coincides 
with the model of Ref.^^ provided the latter is consid- 
ered in the limit of a small angle weak scattering. The 
reduced cr-model obtained in this way is equivalent to 
the model found in the problem of a long range poten- 
tial disorder Ref.^''. At the same time, it is relevant to 
emphasize that at short distances inside the Lyapunov 
region the RMF and RP models correspond to different 
cr-models. 

At distances, exceeding the transport length Itr = 
vpTtr one comes to the standard diffusion cr-model, Eq. 
(3.23), unless the correlation of the magnetic fields obeys 
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Eq. (3.24). Calculations for the tr-modcl, Eq. (3.23), 
within the renormalization group scheme leads to the 
standard conclusion about the localization. This con- 
clusion is in contradiction with the numerical results of 
Ref.^^ where the existence of "hidden degrees of free- 
dom" was proposed, which could lead to the existence 
of extended states. We did not find any indication for 
such degrees of freedom. Of course, our consideration 
was performed in the quasiclassical limit, such that we 
did not take into account a possibility of a quantization of 
the energy levels. However, it is not easy to understand 
how taking into account distances shorter than the wave 
length Xf could lead to a destruction of the localization. 



APPENDIX A: LYAPUNOV EXPONENT IN RMF 
PROBLEM 

Here we study the classical scattering of two particles 
in a random magnetic field (RMF). The presence of the 
RMF leads to an efi^ective interaction between the parti- 
cles. The radius of this interaction is equal to the corre- 
lation length of the field. The scattering process lasts a 
finite time after which the particles diverge over the dis- 
tance exceeding the correlation length and begin to move 
without any interaction. The aim of the calculation pre- 
sented below is to estimate this time. It is clear that for 
larger times the particle scattering may be considered in 
terms of collisions. We restrict our calculation by the 
case of a weak magnetic field. 

Let us consider two particles on a plane with the coor- 
dinates ri, r2 and momenta pi, p2 moving in a perpen- 
dicular magnetic field. The equations of the motion for 
each particle are 



p, . eB{ri) 

— , Pi = — ^ — Pi 
m mc 



X e. 



(Al) 



where is the unit vector perpendicular to the plane of 
the motion. Let p = ri — r2 and p = pi — P2 be coordi- 
nate and momentum of the relative motion. We assume 
that the particles start their motion close to each other 
and have parallel momenta pi = p2 so that p = and 
p = Pq in the beginning; po is assumed to be perpendicu- 
lar to the direction n of the motion of the center of mass. 
Since the energy does not change in the magnetic field, 
the absolute value of the momenta pi , p2 will remain con- 
stant and equal to each other |pi| = |p2|. Therefore, the 
direction of the relative motion will always be perpendic- 
ular to the direction of the motion of the mass center n: 
(pn) = 0. This allows us to write p = p[n x e^]. Using 
Eq.(Al) we find 



p=— , P = e — [B1-B2) 
m c 



(A2) 



where Bi = B{ri), and p = |p| is the absolute value of 
the momentum of the relative motion. At the beginning 
of the motion p is rather small and the diff'erence Bi — B2 



can be approximately written as Bi — B2 ~ (dB /dR±)p, 
where R± is the coordinate of the mass center in the di- 
rection perpendicular to n. Eq.(A2) considered in this 
approximation reduces to a linear system of first order 
differential equations. Hence, the distance p will grow 
exponentially as a function of time. The mean rate of 
the divergency or the Lyapunov exponent determines the 
scattering time involved. 

To study statistics of the relative motion we introduce 
a distribution function W{t,p,p). By definition, it is the 
probability for the relative distance and momentum to 
be p and p at the time t, respectively, provided they have 
been initially po, p = 0- Let W{to,p,p) be the distribu- 
tion at the time to- Then, it can be written at the time 
to + At as 

W{to + At,p,p) = j P{to + At,p,p;to,p',p') 

xW{to,p',p')dp'dp' (A3) 

where P{t, p,p\t' , p' ^p') is the transition probability. 
This probability is determined by the equation of mo- 
tion, Eq. (A2), and is introduced as 

P{t, p,p; t', p',p') = s(^p-p'- £ ^dr^ 

x6(p-p'-e^ [Bi(t) - B2iT)]dT^ (A4) 

where B,{t) = B[r,{T)], r,(T) = R(t)±p(t)/2 and p(t), 
p{t) are the solution of the classical motion equation 
(A2). Substitution of Eq.(A4) into Eq.(A3) gives a rela- 
tion between the distributions W at times to and t^ + At. 
Assiiming that At is smaller than the inverse Lyapunov 
exponent we expand this relation in At and then average 
over the magnetic field B{r). Since the magnetic field is 
assumed to be weak, we neglect the influence of the field 
on the trajectory of the mass center and obtain 



dW 



p dW 
m dp 



2 . ^ 2 d'^W „ 



(A5) 



where Ttr is the transport time, Eq.(3.22), that can also 
be written as 



1 f~^°° 

Ttr J-00 



The function e{p) is by definition 



/_oo WBivpnT + p[n X e^])dT 
J^^ WB{vFnT)dT 



(A6) 



(A7) 



The distance between the particles in the Lyapunov re- 
gion is smaller than the correlation length of the magnetic 
field b. Hence, one may expand the function e{p) in p, 
which gives e(p) w p^/2b^. This relation is to be con- 
sidered as a definition of the length b. Substituting this 
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expansion into Eq.(A5) we come to the same equation as 
the one derived in Ref.^* where electron scattering in a 
long-ranged potential disorder was considered. 



d a p2 

— - VF(P^ 2 

at op Ttrd'^ 
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W = 



(A8) 



Using the result of that paper wo find that the func- 
tion W{t,p) determining the distribution of the distance 
p (the momentum of the relative motion p is implied to 
be averaged in this function) satisfies the equation 



■ d d- 

ot oz 



W = Q 



(A9) 



where /3 is a numerical coefficient equal to /? ~ 0.365 and 
z = \n{b/p). It follows from Eq.(A9) that the coefficient 
tl is in fact a characteristic time of the divergency of the 
trajectories of the particles calculated from the classical 
motion equation Eq.(A2). According to Ref.^^ this time 
is equal to 



ihr) 



2/3 



(AlO) 



and this is at the same time the inverse Lyapunov ex- 
ponent. As mentioned above, the quantity tl has the 
meaning of a characteristic time that two scattered par- 
ticles spend moving together until the distance between 
them starts exceeding the correlation length b. 
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